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Abstract 



We argue that the well-known problem of the instabilities associated 
with the self-forces (radiation reaction forces) in classical electrodynamics 
C3 ' are possibly stabilized by the introduction of 

gravitational forces via general relativity. 



1 Introduction 

The problems and difficulties associated with the motion of charged particles 
interacting with both an external electromagnetic field and its own self-field 
have not been resolved even after over a century of investigation^ [2j O |H [5] . 
The problems arise in many different contexts: the difficulties in giving appro- 
priate initial conditions, infinite self-energy problems, model building, Lorentz 
invariance difficulties and perhaps the most serious, the instabilities (or pre- 
acceleration) in the solutions to the equations of motion. 
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The best known equations of motion, coming from a point structureless 
particle, are the Abraham-Lorentz equations and the relativist generalization 
the Abraham-Lorentz-Dirac equations. They are given, respectively, by[2] 



mv — qE + qB x v + —r v , 




3c ; 



(1) 



and[2] 



mv 



• a 



qF ab v b + ^(v a + ^v a v b v b ), 



(2) 



where F ab (or E,B) are an external held and are derived by a variety of means, 
but always with severe approximations. 

It is generally acknowledged that there are fundamental difficulties with this 
issue. And there seems to be a variety of different reasons, explanations and 
suggested remedies. They range from: quantum theory is the resolution to the 
approximations leading to these equations are wrong or even that there is no 
real problem. The author, J. D. Jackson, summarizes the situation in his well 
known graduate textpQ as: 

"The difficulties presented by this problem touch one of the most funda- 
mental aspects of physics, the nature of the elementary particle. Although 
partial solutions, workable within limited areas, can be given, the basic prob- 
lem remains unsolved. One might hope that the transition from classical to 
quantum-mechanical treatments would remove the difficulties. While there is 
still hope that this may eventually occur, the present quantum-mechanical dis- 
cussions are beset with even more elaborate troubles than the classical ones. It 
is one of the triumphs of comparatively recent years (~1948 - 1950) that the 
concepts of Lorentz covariance and gauge invariance were exploited sufficiently 
cleverly to circumvent these difficulties in quantum electrodynamics and so al- 
low the calculation of very small radiative effects to extremely high precision, 
in full agreement with experiment. From a fundamental point of view, however, 
the difficulties remain." 

The purpose of this note is to describe a totally new point of view towards 
this problem. This new view is completely classical, with no reliance on quan- 
tum theory. It however does rely heavily on the Einstein-Maxwell equations of 
general relativity. 

The basic situation that we address is to hrst consider an arbitrary com- 
pact gravitating-electromagnetic system which is taken to be the particle whose 
motion we want to describe. The system is given by local mass and charge 
densities and currents. There are no external helds acting on it. It is an iso- 
lated system with arbitrary internal degrees of fredom. The program is to solve 
the Einstein-Maxwell equations in the future null asymptotic region and, from 
the asymptotic field (the asymptotic Weyl and Maxwell tensors), determine a 
center of mass and center of charge and their laws of motion. Aside from condi- 
tions that the total charge Q be non-vanishing and the important requirement 
that the centers of mass and charge should coincide, there is no further model 
building. 
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Since our detailed calculations, which were done in the language of the spin- 
coefficient formalism, are long and complicated and have appeared elsewhere [6], 
we will just summarize the ideas and results. Basically the calculations are done 
2 nd order in deviations from Reissner-Nordstrom. In the spherical harmonic 
expansions, with frequent use of Clebsch-Gordon products, only terms up to 
the I = 2 harmonics are kept. 



2 The Complex Center of charge: Its Identifica- 
tion 

The basic starting idea in this work is essentially simple. It is in the generaliza- 
tions and implementations where difficulties arise. 

Starting in Minkowski space in a given Lorentzian frame with spatial origin, 
the electric dipole moment D e is calculated from an integral over the (localized) 
charge distribution. If there is a shift, R, in the origin, the dipole transforms 
as 

Tr E = D E - QR. (3) 

If D e is time dependent, we obtain the center of charge world- line by taking 
D* E = 0, i.e., from R = De/Q- It is this idea that we want to generalize and 
extend to gravitational fields. 

First, however, we want to discuss other dipole issues in flat-space. Starting 
on the time-like world-line at the spatial origin, we construct the family of future 
directed light-cones, Co, and investigate behavior of the Maxwell field in the limit 
as null infinity is approached, i.e., Penrose's 3 + . Using the null tetrad formalism 
and where the Maxwell pair (E , B or F a i,) is replaced by the complex vector 
E + i-B, or more accurately by their tetrad components, (0o, 01, 4>2), with[7] 

00 = F ab l a m b (4) 

01 = ^F ab (l a n b + m a m b ) 
4>2 = F a bffi a n b . 

The asymptotic (peeling) behavior of these fields for a compact source is given 

by 

^ = $ + °(^ 4 ) (5) 
- - ^ ■ 0(r -») 



•J2B 



r 



0(r- 2 ). 



The vector field l a is the tangent field to the null geodesic generators of the 
null cones Co- At 3 + , n b is the tangent field to the null generators of 3 + while 
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(m a ,m b ) are (the complex conjugate pair) tangent to the two surface, S 2 , the 
intersection of Co with J + . 

The r independent quantities (0q, 0?> $2' •••) are functions 'living' on 3 + , 
i.e., functions of the retarded time, u, (the light-cone cuts of 3+) and (C,C) ; 
the complex stereographic coordinates labeling the, S 2 , generators of 3 + . The 
components of their spherical harmonic decomposition, 

K = 4iY 1 \ + 4>°o ij Y 2 1 ij + ..., (6) 
$ = Q + *M + + (7) 

^^^ t Yu 1 +^Y^ + ..., (8) 

are the asymptotically defined multipole moments and their time derivatives. 

For example, the I — harmonic component of is proportional to the total 
source charge. For us the important quantity is cj)^, the I — 1 component of 
4>q '■ is proportional to the {asymptotically defined) complex dipole moment, 

D c — De + iDM, where Dm is the magnetic dipole moment. 

The problem now is: how does the Dq transform under an origin shift to 
an arbitrary world-line? With an origin shift there will be new light-cones and 
a new null vector field, l* a , obtained from the old one, l a , by a null rotation at 
J + . This can be expressed explicitly bv[9lfT0] 

I* = I + -m + -m + 0(r- 2 ) (9) 
r r 

777* = 771 + 0(f" ) 
77 = 77 

where L = L(u, (, £) is a stereographic angle field given on 3 + (still to be 
described) that determines the new null geodesic field, I*. 

The transformation law for <j}^ i (given only approximately for small origin 
shifts [8]) is® 

<roi — wo 

If we are given a Minkowski space world- line, x a = £, a (s), for the apex of the 
new light cones, then L = L(u, £, Q is given in the parametric form 

L(u,c,c) = r(s)TO (c,c), (ii) 



*°-^° 2LcjP x +,..) i (10) 



with 



ycO^dM,-^,^), ( i2) 
y ^ ' 2 v l + CC i + CC i + CC 

m a (C,C) = (O.l&CC)) = ^(0,l-f,-7(l + C 2 ), 2C), (13) 



By the appropriate choice of £ a (s), from Eq. ([i"0|) , with the use of Ea. (fTT|) . 

y o i 



one can force the real part of (j)^ to vanish, thereby making x a = £ a (s) the 
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center of charge. If however we generalized the choice of L(u, £, C) and allowed 
it to be defined parametrically by 



£(«,C,0=&(r)m o (C,0. ( 14 ) 

where z a = is a complex analytic world in complex Minkowski space, then by 
setting = 0, in Eq.flOf the complex world-line is determined. This complex 
curve (which is purely formal) defines the complex center of charge. Using this 
"curve" as the origin both the elecric and magnetic dipoles vanish. 
Theorem [11 J - copy from home 



3 The Complex Center of Mass 

For asymptotically flat Einstein-Maxwell space-times the situation is totally 
analogous: the shear-free null geodesies originating from light-cones from world- 
lines (real or complex) are replaced by (regular) asymptotically shear-free null 
geodesic congruences generated by a complex world-line 11 , in the space of the 
complex Poincare translation subgroup of the BMS group. The Maxwell asymp- 
totic dipole transforms exactly as in the flat space case, i.e., as in Eq. (fTTj|) with 
however a slight change in the parametric description of the function L(u, C, C) : 



L = C{r)Y^CX) - 6e j (r)Y 2 \j(C,0, (15) 
u = -±=e(r) - \C (r)r°(C,C) + C'(r)Y 2 %(CX) + (16) 

where the extra terms come from the existence of a non-vanishing Bondi shear, 
given up to I = 2 terms, by 

a = 2AC\u)Y 2 % + ... (17) 

Turning the gravitational behavior, the relevant (for us) tetrad components 
of the Weyl tensor [TUl [5] 



V'i = -C abcd l a m b l c m 
tp 2 = -C a bcdm a n b l c rr 

have the asymptotic form (the peeling theorem) 



V 
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The leading terms have the harmonic expansion: 



^2° = T + + + ... (18) 

^i=^+<-^ + - (19) 
The mass aspect, defined by 

is rea/ and has the expansion 

* = * = *° + + ¥ j Y^ + ... 

The Bondi mass and linear momentum (four-momentum) is obtained from the 
I = (0, 1) harmonic components of W by: 

*° = - 2 -^M (20) 
<r 

^i = _^Lpi (21) 
c A 

The complex gravitational dipole moment (roughly mass dipole plus i times an- 
gular momentum) is identified as being proportional to the I — 1 harmonic of 
tpi, i.e., as (Though many authors add further terms to tp^ for this identi- 
fication, to our order of approximation they all agree with our identification [6].) 

The transformation (to second order) of tp^ to an arbitrary (complex) world- 
line, analogous to Eq. (fTU|) . using Eq. ([T5)) . is[Sj 

< = + ...)«• (22) 

Setting ipl® = 0, thereby defining the complex center or mass, yields after 

a lengthy calculation, 

= _^2G M[rH + il ekji v k ¥] + G\ (23) 

where G % is a known non-linear function of quadrupole terms. 

At this point we make our only assumption on the physical system being 
considered. We saw that we could determine a complex center of charge or a 
complex gravitational center of mass by setting either ip^ or ip^ to zero. We 
now assume, for the rest of this work, that the two complex world-lines coincide. 
Aside from taking Q ^ 0, there are no other conditions on the internal structure 
of our source (particle). 

We now turn to the dynamics, which are contained in the asymptotic Bianchi 
identities. They can be written: 



G 



(V>? )■ = -m + Scr(ef)- + 3ad{wy + d 3 W + 2kfy 2 (24) 

= aa' + fc^202 ( 25 ) 
k = 2Gc~ 4 (26) 



From these two equations, (f2"4")l and (|23|) . we extract the equations of motion 
with the radiation reaction term. Rather than going thru the details (long with 
rather unattractive calculations) we will describe what we did in words and 
then give the results. First we point out that we change the Bondi -^-variable 
to w = -\Z2uc~ 1 , w being the retarded time. Derivatives with respect to w are 
denoted by prime, ('). 

We first extract from Eq. (|2"4]) its / = 1 part and then decompose it into its 
real and imaginary parts. This yields two results: the imaginary part determines 
the dynamics of the total angular momentum, i.e., the conservation of angular 
momentum. Other than remarking that we identify S l — Met;} as the intrinsic 
spin (with the imaginary part of £ l ), this is not our interest here and will 
not be discussed any further. The real I = 1 part can be solved for the linear 
momentum P l that was sitting in the I = 1 part of : 

P k = Mv k R - ||-t# + W k . (27) 

This is a major result that come from our identification of the complex centers 
of mass and charge. First of all we see kinematic expressions for the Bondi 
3-momentum, the mv term and then the radiation reaction contribution to 
the momentum. The W contains further kinematic terms involving spin and 
quadrupole interactions that are known but not displayed here[B]. 

Extracting the I — (0, 1) harmonics from Eq. (l2"5l) yields the Bondi mass and 
momentum loss equations: 

M ' = -^m'LQZ'L + Q'iXL) - + <<) (28) 

180c 7 1 

(29) 

4. (J ' n jk "' + v 3 ' D jk "') + 1 ( D lj "' D ij "' - D lj "' D ij "')f,i, 
15c 5 M > 54Q C 6 ^ E M M e 

with the mass and spin quadrupoles related to the ^ by 



e j = + «#) - y^(QmL s + iQ%J. (30) 
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The mass loss equation is thus exactly the usual quadrupole energy loss plus 
the classical dipole and quadrupole electromagnetic energy loss. 

It is however Eq. (|2T))) that is of most interest to us. By substituting the 
kinematic expression for the momentum, Eq. (|27p into Eq. (|29p we obtain our 
generalized Abraham-Lorentz equations of motion: 

Mv% + v k R M' -^v k R " +R k =F k . (31) 

Note that though it is similar to the Abraham-Lorentz equations there are 
many differences that are hidden in the known but complicated expressions for 
M', R k and F k . The F k is the Bondi recoil (or rocket) force due to the momen- 
tum loss, while R k can be considered to be a gravitational radiation reaction 
force depending on internal degrees of freedom spin and quadrupole moments. 
M' has exactly the classical dipole energy loss terms plus two additional terms. 
Though it is very hard to directly see if the solutions to Eq. (|3"Tj) are well behaved, 
in the conclusion we will discuss this issue in more general terms. 



4 Conclusions 

We have considered the situation of a gravitating - electromagnetic source of 
compact support viewed from future null infinity. The only restriction made on 
the distributions is that the total charge is non-vanishing and that the complex 
asymptotic electromagnetic dipole moment be proportional the complex gravita- 
tional dipole moment so that the complex centers of charge and mass coincide. 
Though it is not clear how severe this condition is, it certainly is a serious re- 
striction. It has been shown that for this situation the gyromagnetic ration, the 
ratio of the spin-angular momentum to the magnetic moment is that of Dirac's, 
namely g — 2. We showed that in a manner completely analogous to the flat 
space Maxwell case, one could determine the transformation laws for the two 
dipole moments and thereby go to the center of mass/charge, determining a 
unique complex world-line. Then, using the Bianchi Identities, which play the 
role of dynamical equations, we were able to give kinematic significance to the 
Bondi linear momentum, in the sense that we had 

P = Ml? -^rV + ... (32) 
3c 6 

From the Bondi momentum loss equation it immediately followed that we 
had a generalized version of the Abraham-Lorentz equations of motion for an 
isolated massive charged particle. It should be emphasized that the quadrupole 
quantity, £ lJ (u), is arbitrary and in most cases it is taken as non- vanishing in a 
finite interval so that the gravitational radiation exists also in a finite interval. 
If however the motion, from Eq. (|31[) . is unstable, the particle acceleration will 
be unbounded and there will be an infinite energy loss via the electromagnetic 
dipole radiation. This would be a physically unacceptable situation, indicating 
that something is seriously wrong with the Einstein-Maxwell equations. 
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The question then is does the general relativity (gravitational) contribu- 
tions to the equations of motion stabilize the equations. Though we do not see 
any immediate prospects for a direct proof, we make a few comments. Looking 
at Eq. (|3"Tj) . we see that the term M'v has the same form as in the Abraham- 
Lorentz equation but now is more negative because of the extra radiation terms 
and has the correct sign to try to stabilize the motion. Whether or not it does 
stabilize is an open question. If not, perhaps higher order terms that have been 
left out in our approximations might succeed. And finally, it is known that the 
vacuum Einstein equations are stable in the neighborhood of Minkowski space. 
If the same were true of the Einstein-Maxwell equations with compact sources, 
that would constitute a proof that our physical situation was indeed stable and 
the run-away behavior was prevented by the inclusion of classical general rela- 
tivity. Unfortunately, this is a difficult question and, to our understanding, the 
answer is unknown. It would be surprising if it turned out that the asymptotic 
Einstein-Maxwell equations were unstable. 
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